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Ni 8 oFe 2 o (Py) and Py-Cu exhibit intriguing ultrafast demagnetization behavior, 
where the Ni magnetic moment shows a delayed response relative to the Fe [S. 
Mathias et al., PNAS 109 , 4792 (2012)]. To unravel the mechanism responsible for 
this behavior, we have studied Py-Cu alloys for a wide range of Cu concentrations 
using X-ray magnetic circular dichroism (XMCD). The magnetic moments of Fe 
and Ni are found to respond very differently to Cu alloying: Fe becomes a strong 
ferromagnet in Py, with the magnetic moment largely unaffected by Cu alloying. 

In contrast, the Ni magnetic moment decreases continuously with increasing Cu 
concentration. Our results are corroborated by ab-initio calculations of the electronic 
structure, which we discuss in the framework of virtual bound states (VBSs). For 
high Cu concentrations, Ni exhibits VBSs below the Fermi level, which are likely 
responsible for an increased orbital/spin magnetic ratio at high Cu concentrations. 

Fe exhibits VBSs in the minority band, approximately 1 eV above the Fermi level in 
pure Py, that move closer to the Fermi level upon Cu alloying. A strong interaction 
between the VBSs and excited electrons above the Fermi level enhances the formation 
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of localized magnons at Fe sites, which explains the different behavior between Fe 
and Ni during ultrafast demagnetization. 
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I. INTRODUCTION 

Ni 8 oFe 2 o (Py) is a commonly used material for studying effects involving spin polarized 
currents, such as spin transfer torque |T], magnetoresistance [1]-|3], spin pumping|lj, spin hall 
effect [3], inverse spin hall effect [5|, and ferromagnetic Josephson junctions!?]. It is not only 
used for its low magnetostriction (no strain induced by magnetization) and the combination 
of high permeability and low coercivity, but also for its high spin polarization of conduction 
electrons and large resistivity difference between minority and majority spins[3l |8|. The 
origin of this spin dependent resistivity comes mainly from three different effects: magnon 
scattering, phonon scattering and impurity scattering [9]. The spin dependent impurity 
scattering, which has a major contribution even at room temperature, comes from the large 
impurity potential of Fe in the minority band of Ni[2l ITU] . 

In a rigid band model (RBM) [TT] , where the constituents of an alloy are assumed to form a 
shared band structure, these types of impurity potentials would not exist. It has been shown 
both experimentally [12] and theoretically that the RBM is inaccurate for transition metal 
(TM) alloys, but the behavior of these alloys is anyway described in the spirit of the RBM in 
many textbooks [T3ll2^ . Consequently, the RBM has often been used to explain the Slater- 
Pauling (SP) curve, while the formation of virtual bound states (VBSs) was successfully 
used for explaining its deviating branches that exhibit antiferromagnetic interactions, such 
as Cr and V alloyed with Co and Nij^S] [23]. The formation of VBSs for dilute alloying 
concentrations has been found theoretically for a wide range of TM alloys [2?ll37j . Although, 
ab-initio calculations have been performed for Py in several publications [321 jSSj-HO], there 
are basically no discussions about the most prominent feature in its density of states (DOS) 
and its implications, i.e. the formation of virtual bound states for Fe minority d-states, 
which form due to the relatively low concentration of Fe (20%) and the energy mismatch 
between the minority d-states of Ni and Fe. 

We present x-ray magnetic circular dichroism (XMCD) measurements on a series of 
(Nio, 8 Feo. 2 )i-xCu 3 , (i.e. Pyi_ 3 ;Cuj;) alloys, where we study the charge transfer and the ele¬ 
ment specihc magnetic moments of Fe and Ni. We show the existence of VBSs in Pyi_ 3 ;Cu 3 ; 
by XMCD, and study their behavior using ab-initio calculations. We use these results to 
shed light on a currently unexplained ultrafast phenomena in Py, where the magnetizations 
of Fe and Ni show very different responses to excitations by an intense laser pulse. Mathias 
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et al.[lT], found that the Ni magnetization appears to be initially unaffected by the laser 
pump light and instead follows the demagnetization of the Fe, which responds quickly to 
the laser light, with an 18 fs delay in Py and 76 fs delay in Py 6 oCu 4 o. Considering that pure 
(not alloyed) Ni and Fe demagnetize on similar time scales, their conclusion might seem 
surprising. Nevertheless, the results have been reproduced by others [4^ 148] and alternative 
explanations have not emerged. Importantly, a common misconception is that results by 
Eschenlohr [44] contradict these Endings. However, as explicitly explained in Ref. |44], a large 
uncertainty in the determination of the pump-pulse arrival time was handled by shifting the 
Fe and Ni data sets to eliminate any time-delay between Fe and Ni. 

Fe is found to be a strong ferromagnet in Py, i.e. the Fe majority d-states are filled when 
disregarding s — p — d hybridized states above the Fermi level, and its magnetic moment is 
largely unaffected by Cu alloying. Ni shows a continuously decreasing magnetic moment with 
Cu alloying. Ab-initio calculations support these observations and also show the formation 
of VBSs in the Py DOS. Ni has VBS’s below the Fermi level for high Cu concentrations and 
Fe exhibits VBSs in the minority band, ~1 eV above the Fermi level in pure Py, which move 
closer to the Fermi level with Cu alloying. The resulting strong interaction between the 
VBSs and conduction electrons enhances the exchange interaction between Fe sites. Most 
interestingly, our results support a picture where Elliot-Yafet spin-flip scattering |441 SS] is 
only indirectly responsible for the demagnetization in Py, through an ultrafast generation of 
magnons. The magnon generation is highly promoted at the Fe sites due to the formation 
of VBSs and hence the demagnetization of Ni depends on the time it takes for the magnons, 
which are initially localized at Fe atoms, to distribute throughout the Py lattice. This 
mechanism for demagnetization is fully consistent with the results and conclusions given by 
Mathias et al.|41j. 


II. EXPERIMENTAL AND THEORETICAL METHODOLOGY 

Thin films of 20 nm (Feo. 2 Nio, 8 )i-xCua;, with 0 < x < 0.55 were grown on 3 nm Ta seed 
layers by magneton sputtering on Si 3 N 4 substrates. The samples were capped by 3 nm Ta 
to prevent oxidation. Two sample sets, SI and S2, were grown at two different occasions in 
the same sputtering system. The results from both sample sets were practically identical, 
confirming the repeatability of our measurements. Here we present the averaged data of 
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the two samples sets. These samples were grown in the same spattering system and with 
the same growth parameters as for the samples stndied by Mathias et al.|lT] (see their 
supporting information for sample characterization). The XMCD measurements were per¬ 
formed at beamline U4B at the National Synchrotron Light Source at Brookhaven National 
Laboratory in (Upton, NY), using 70% circularly polarized light at normal incidence to the 
samples. The samples were saturated out-of-plane (OOP) using a magnetic held of 1.5 T. 
The magnetic contrast was obtained by switching the OOP magnetic held, where the cor¬ 
responding spectra are referred to as M-|- and M—. The measurements were performed in 
transmission. Unlike total electron yield (TEY) measurements, which are more commonly 
used, the transmission geometry is insensitive to externally applied helds and is not ahected 
by saturation ehects which can ahect the derived magnetic moments 071 . Also, transmission 
measurements provide a much more reliable method for comparing the absorption coefficient 
between diherent samples, as they are insensitive to sample dependent electron escape depth 
and generation of secondary electrons [17]. A reference sample without the Cu-Py layer was 
used to correct for the absorption from the Si 3 N 4 and Ta layers. The provided uncertain¬ 
ties correspond to the standard deviation and are derived by estimating the variance in the 
diherence and sum signal of the M-|- and M— spectra. Also, an estimated uncertainty in 
the htting of the step function was included in the total variance. The accuracy of the him 
thickness and Cu concentration will affect the analysis of the total absorption cross-section, 
and was estimated to have a standard deviation of 3%, as determined from x-ray reflectivity 
measurements. 


Ab-initio calculations were performed using the spin-polarized relativistic Korringa-Kohn- 
Rostoker method in combination with the coherent potential approximation fCPA) |l8l 
To investigate the effect of the exchange-correlation functionals, we used several different 
implementations of the local density approximation and generalized gradient approximation. 
We found no signihcant dependence on the exchange-correlation functional for the quantita¬ 
tive values. Hence, we only present the values obtained from a full potential spin-polarized 
scalar relativistic calculation with the Perdew-Burke-Ernzerhof functional [50 j for exchange 
correlations. 
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FIG. 1. Spin magnetic moment of Ni and Fe, respectively, as a function of Cu concentration. 
The theoretical (black dashed) Ni moment decreases linearly with Cu concentration while the Fe 
moment remains constant. Inset in a) shows the absorption and XMCD of Ni in Py to illustrate 
the parameters used in Eq. Inset in b) shows the nearest neighbor Fe-Fe exchange interaction 
from ab-initio calculations, and an estimated RKKY-like exchange between magnetic 

impurities 


III. EXPERIMENTAL RESULTS 


The spin and orbital magnetic moments are given by the XMCD snm rnles[5l] and are 
proportional to the number of 3d-holes {Nh), 

AL3 — 2AL2 

- ■ 

( 1 ) 


m., = 


2 AL3 + AL2 

mi = - - - Nh 

6 r 

The parameter r is the total integrated intensity of the absorption edge (averaged over both 
magnetization directions) after subtracting a background (see Fig[^, inset). The position of 
the steps in the background were taken at the energies corresponding to the peak absorptions 
of the L 3 and L 2 edges with a step height ratio of 2:1. AL„ (n = 2,3) are the integrated 
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absorption differences between M+ and M- spectra at the absorption edges, as illustrated 
in Fig. inset. The magnetic dipole operator can safely be ignored since the sample is thick 
enough to make any interface effects negligible |S2]- The spin moments obtained from XMCD 
and ab-initio calculations are presented in Fig. [T^ and b for Ni and Fe, respectively, for 
different Cu concentrations. The ab-initio spin moments are shown as black dashed/circles, 
which correspond to a 0 K magnetization. The green dashed/triangles correspond to ab- 
initio values after correcting for a decreased magnetization at room temperature. The room 
temperature (RT) correction was obtained from SQUID magnetometry, where we measured 
the ratio in total magnetization between 300 K and 10 K. Note that this correction is fully 
valid only when Fe and Ni have equivalent temperature dependence. Since both Fe and Ni 
show an almost identical XMCD ratio between 300 K and 80 K for 40% Cu, we can assume 
this approximation to be reasonable. For deriving the XMCD spin moments, we used the 
number of holes in Eq. as a htting parameter to equate the spin magnetic moment from 
XMCD and ab-initio at RT for permalloy (Py, 0% Cu), giving 3.41 holes for Fe and 1.28 
holes for Ni in Py. 

The spin moments in Py correspond well to values found by neutron scattering techniques 
|53j . where the magnetic moment for Ni is similar to the pure element (~ 0.6/rs) and the 
magnetic moment of Fe is enhanced to about ~ 2.7fis compared to ~ 2.2/is in pure bcc Fe. 
The theoretical Ni spin moment in Fig. [pi shows an almost linear decrease with increasing Cu 
concentration but remains non-zero even up to 90% Cu. The RT XMCD Ni moment shows 
values that are somewhat higher than what is obtained from RT-corrected ab-initio values. 
However, the Ni moment follows a similar trend with Cu concentration. It is possible that 
the number of Ni 3d-holes used for deriving the magnetic moment has been overestimated, 
since a lower number of holes would provide a better correspondence between theory and 
experiment for Cu alloyed samples. There have been speculations about the mechanism 
for the decreasing Ni moment in Nii_ 3 ;Cua: alloys; however, no clear understanding has yet 
emerged [SIHSS]- Our conclusion is, with the support from the calculated DOS discussed 
later, that the localization and interaction between Cu and Ni d-states drives the reduction 
of the Ni magnetic moment. For the Fe spin moment shown in Fig. [p), theory suggest an 
almost constant value that persist all the way up to 100% Cu at T = 0 K. The RT XMCD 
measurements follows closely the RT-corrected ab-initio values, strongly supporting the 0 
K ab-initio result that the local Fe magnetic moment is almost unaffected by Cu alloying. 
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FIG. 2. The ratios between the orbital and spin moments for Fe and Ni in Pyi_a;Cua; as a function 
of Cu concentration. The ratio remains constant for low concentrations of Cu alloying. At high 
Cu concentrations (> 40%) the ratio increases for Ni. 

The 80 K magnetization is slightly lower than expected from SQUID measurements. This 
could possibly be attributed to a higher sample temperature than the temperature reading 
on the sample holder. 

The orbital moment has a very similar dependence on Cu concentration as the spin mag¬ 
netic moment. To identify any Cu concentration dependent trends in the orbital moment, 
we studied the orbital/spin ratio, shown in Fig. Unlike the spin and orbital moment by 
themselves, the mi/rris ratio is independent of the number of Sd-holes, as apparent from 
Eq. The values are slightly higher than what has been measured in other studies for bcc 
Fe (0.043)[5l], Ni (0.11)[57] and Py (Ni; 0.10, Fe: 0.04)|5H]. The ratio for both Fe and Ni 
remains almost constant through a wide range of Cu concentrations. However, there appears 
to be an increase in the ratio for Ni at high Cu concentrations. This is likely related to the 
narrowing of occupied Ni d-states, see inset in Fig. This will be discussed later in relation 
to the calculated density of states. 

The intensity of the L 2 3 absorption edges corresponds to electron transitions from 2pi 3 

’ 2 ’ 2 

core levels to empty 3d states and is proportional to the number of 3d-holes|59]. In Fig. 

and b, we show the variation in the number of 3d-holes obtained from the absorption 
intensity of Ni and Fe, respectively, as a function of Cu concentration. We choose to show 
both sample sets (SI and S2) separately since, unlike XMCD, these values are sensitive to 
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FIG. 3. The number of Sd-holes in a) Ni and b) Fe, obtained from the white-line intensity of the 
L 2,3 absorption edges. The number of Fe 2>d holes appears to be unaffected by Cu alloying while 
there is a small charge transfer to Ni Sd-states. 


variation in both sample thickness and elemental concentrations. The numbers of Sd-holes 
obtained from our calculations are shown as black dashed lines. The number of holes in 
Ni, shown in Fig. obtained from ab-initio is about 10% higher than what we obtain 
from the sum rules, but both experiment and theory show a very similar decreasing trend 
in the number of holes with increased Cu concentration. This indicates either a charge 
transfer from Cu or a charge reordering between s — p and d-states in Ni. In either case, 
the decrease is very small and can only account for a decreased Ni magnetic moment of 
10% (for 55% Cu) at most. Similar values have been found experimentally for Nii_ 3 ;Cu 3 ; 
alloys[Hn] (~ 5% decrease of Ni d-holes with 55% Cu). Theoretically, Stefanou et al.[M] 
has provided a detailed calculation of charge effects for impurities in Ni, which corroborates 
the weak charge transfer between Cu and Ni. The number of holes in Fe, shown in Fig. 
i). gives a much smaller deviation between the XCMD and ab-initio values than Ni. Both 
theory and experiment indicate an almost constant value of Fe 3d-holes, independent of the 
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Cu concentration. 

The high magnetic moment of Fe, its insensitivity to Cu alloying, and the difference in 3d- 
state Filing between Fe and Ni, lead us to the following conclusions about the Fe impurities: 
(1) Since Fe has approximately 2 /i^ higher magnetic moment than the strong ferromagnet 
Ni, which corresponds to the difference in d-state filling (2e“, see Fig. between Fe and 
Ni, then only the filling of the minority band can differ between Fe and Ni in Py. Hence, 
Fe behaves as a strong ferromagnet, with a filled majority band in Py, regardless of the Cu 
concentration. (2) Since Ni and Cu have nearly filled d-bands, the relatively large number 
of empty minority d-states in Fe should not hybridize well with the Ni and Cu d-bands. (3) 
The exchange splitting of the Fe atoms is largely insensitive to the exchange splitting of the 
Ni-host, as evidenced by the insensitivity of the Fe moment to Cu concentration in contrast 
to the Ni behavior. From the above statements one can conclude that the minority d-states 
of Fe are decoupled from the host minority d-band and could hence be described as virtual 
bound states above the top of the host d-band, as will be described in the following section. 


IV. VIRTUAL BOUND STATES (VBS) 

The valence band of 3d-metals consists of s, p, and d-states. Electrons in s and p states 
are loosely bound to atoms and can easily form electronic bands, resulting in electrons that 
are highly mobile. The d-states can be considered more localized and less mobile, but their 
wave functions still extend far enough to hybridize with d-states of neighboring atoms and 
form d-bands. However, band formation not only requires a spatial overlap of the wave 
functions, but also an overlap in momentum space, i.e. the closer electronic states are in 
both momentum and energy, the stronger they can hybridize. For 3d impurities in a 3d metal 
host, the free electronic s, p-states of the host and impurity will easily form a common s — p 
band. However, the d-states will experience a potential change at the impurity atom, where 
the potential difference between a host and impurity atom is called the impurity potential. 
As illustrated in Fig. |^, the impurity potentials for minority and majority bands, and 
Vf respectively, depend on both Coulomb and exchange interactions |62]. The Coulomb 
contribution is the same for both spin bands. However, the exchange interaction depends on 
the number of occupied states in each band and hence the impurity potentials and I 4 can 
differ. Since the total occupation number is determined by charge neutrality at the impurity 
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FIG. 4. Illustrations of impurity potentials and virtual bound states. (A) The spin dependent 
impurity potential is non-zero at impurity sites and can differ between majority (F]-) and minority 
bands (I^). The impurity potential depends on both the Coulomb and exchange interaction, where 
the Coulomb contribution is equal for both spin bands while the exchange contribution depends on 
the occupation number of each spin band. (B) Without hybridization to s,p-states, the impurity 
potential can induce bound states (BS) at the impurity site. (C) The bound states hybridizes with 
s,p-states and form virtual bound states (VBS). This hybridization induces a broadening of the 
bound states. 

site, the impurity potential Vi will decrease if exchange favors a higher Vf potential, and vice 
versa. For a weak impurity potential, hybridization will occur within the d-band of the host. 
For repulsive impurity potentials i.e. impurities with fewer valence electrons than the host, 
the density of states for the 3d-band will increase at higher energies and decrease at lower 
energies[63]. If the repulsive impurity potential is strong enough, it can even move d-states 
to energies above the d-band, at the impurity site, as illustrated in Fig. |^. These states will 
adopt the character of localized bound states (BS), due to lack of overlap with the host d- 
band [6^l6i] . However, the s—p bands span a very broad energy range. Thus, there is always 
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an energy overlap between 3(i-states and free-electron (s, p) states. The localized Sd-states 
of an impnrity will hybridize with free electronic states and can no-longer be considered 
as bonnd states, bnt are referred to as ’’virtnal bonnd states” (VBS), as illnstrated in Fig. 
|^[65]. Bonnd states have a well defined energy bnt the interaction with free-electronic states 
provide a Lorentzian energy broadening of the VBSs, indicating the non-inhnite lifetime of 
these states. Dne to the interaction between the impnrity 3d-states and free-electron states, 
the VBSs will act as strong scattering centers for condnction electrons if they are close to 
the Fermi level |66j. 


V. SPIN-DEPENDENT DENSITY OF STATES 

In Fig. we present the calcnlated minority and majority 3d-band DOS for a) fcc-Ni, 
b) fcc-Fe, c,e,g) Ni in Pyi_ 2 ,Cna;, and d,f,h) Fe in Pyi_a,Cn 2 , for a variety of different Cn 
concentrations. Both fcc-Ni and fcc-Fe are fonnd to be strong ferromagnets with almost hlled 
majority bands (Fermi level is indicated by vertical dotted lines). The exchange splitting is 
0.7 eV for Ni and 2.5 eV for Fe. Also, the majority and minority bands do not show any 
signihcant qnantitative differences in strnctnre. In Py, the majority states for both Fe and 
Ni are similar to the pnre elements, bnt the minority band of Fe is significantly different, 
see Fig. it- The low energy section (-4 to -6 eV) of the Fe DOS now has a strnctnre 
and exchange splitting which is similar to Ni as observed by comparing the majority and 
magnified (x3) minority band in Fig. [^. This is becanse Fe acts as an impnrity in the 
Ni host band strnctnre with impnrity potentials, Vf and Vi, for the majority and minority 
bands, respectively. As discussed in Appendix B, the main part of the impurity potential 
goes to the minority band for Fe impurities in Ni. This is clearly observed in the minority 
band of Fe, where a large VBS formation is observed, while the disturbance in the majority 
band is insignificant. 

With 20% Cu in Py, both the width and exchange splitting decreases for majority and 
minority states in Ni, as shown in Fig. [^. Since Cu hybridizes with the Ni host and forms a 
joint 3(i-band (see Cu 3d in Fig. |^), it will push Ni 3d states to the higher energy side of the 
3(i-band while Cu mainly occupies the lower energy section. This type of distribution is a 
general feature when two states with different energies hybridize |67j. the lower energy states 
(Cu in this case) will mainly contribute to the lower energy hybridized states, and vice versa. 
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FIG. 5. Spin resolved d-DOS of a) fcc-Ni, b) fcc-Fe, c,e,g) Ni in Pyi-ajCux, and d,f,h) Fe in 
Pyi-xCua; for different Cu concentrations. 

This will also effectively decrease the width of the Ni Sd-DOS and increase the localization 
of Ni 3d electrons relative to pnre Py. For higher Cn concentrations (60%), shown in Fig. [^, 
the localization is so prominent that we observe a well developed VBS below the Fermi level 
in both majority and minority states of Ni. This localization can affect the spin/orbit ratio 
in varions ways and explain the increase observed for Ni in Fig. The spin-orbit conpling 
parameter (^) increases with decreased orbital size and increased gradient of the atomic 
potential, which are both possible conseqnences of VBSs. Also, the crystal field splitting, 
which can originate from either electrostatic interaction between atoms or from covalent 
interactions, should decrease due to the diminished hybridization for VBSs[67l |68]. Any 
combination of these effects will increase the orbital magnetic moment since in perturbation 
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theory the spin-orbit interaction enters the wave-function effectively as (£)"■, where n is 
the order of the perturbation theory. In this expression A is the splitting of the energy 
levels, e.g. due to the crystal held splitting [69]. These effects should not be prominent in Fe 
because all the occupied states are non-localized and only the unoccupied minority states 
show VBS behavior. When the Cu concentration increases, the peak in the DOS for the 
Fe VBS increases, see Fig. |^. This is due to an even greater difference in the number of 
valence electrons between the host and the Fe impurity that results in an increased 1^. The 
black solid vertical line that cuts through Figs. |^, f, and h, shows the energy shift of the 
VBS peak position. The peak position should move to higher energies with increased Vi- 
However, this effect is compensated by a shift of the host minority d-band to lower energies 
due to high Cu content and a reduced Ni exchange splitting. The end result is that the Fe 
minority VBS peak shifts 0.4 eV lower in energies as the Cu concentration increases from 
0 % to 60%. 

The exchange splitting for Ni 3d-states vanishes at high Cu concentrations while it re¬ 
mains large for Fe. This effect was explained by Anderson [65], who showed that the frac¬ 
tional hlling of the d-states together with the Coulomb and exchange interaction determines 
if the ground state is magnetic. In Appendix A, we use the Anderson impurity model to 
demonstrate that Ni and Fe should indeed be non-magnetic and magnetic, respectively, in 
a Cu host. 


VI. FREE ELECTRON-VBS INTERACTION 

An important consequence of the VBS formation is the strong interaction with free elec¬ 
tron states. In the inset of Fig. [^, we show the dependence of the exchange interaction 
upon Cu concentration between two Fe atom when they are nearest neighbors, as 

obtained from our ab-initio calculations. Unlike all other exchange interactions between 
Fe-Ni and Ni-Ni that decreases with Cu concentration, increases strongly. This is a 

direct result of VBS formation in the Fe minority band and their ability to be strong scat¬ 
tering centers for electrons in the minority s — p band. This scattering results in a strong 
exchange interaction between the localized Fe 3d-magnetic moments and the free electrons. 
As we have shown in the previous section, Cu alloying increases the density of VBSs and 
moves them closer to the Fermi level. Therefore, the scattering of conduction electrons at 
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the Fermi level and their exchange interaction with Fe 3d states also increases with Cn con¬ 
centration. This RKKY-like exchange will conple Fe atoms with each other as evidenced in 
the increasing nearest neighbor exchange interaction which is shown in the inset of 

Fig.[p). In the inset, we also show how the RKKY-like pair exchange interaction is expected 
to depend on the Cn concentration in the Anderson impnrity model {J^-k ~ )i where 

ed is the peak energy of the VBS and eu is the energy of the free electrons [701 - 172] . This model 
is described fnrther in Appendix B. Note that this model assnmes highly dilnte localized 
moments in a non-magnetic host. As snch, we only show it here to illustrate that the general 
trend of is almost identical to what is expected from an RKKY-like exchange. Ex¬ 

perimentally, it has been shown that FeCu alloys can indeed exhibit spin-glass behavior |73]. 
which is indicative of such an RKKY-like exchange interactions |67j. 


VII. DISCUSSION 

We have discussed impurity potentials in the framework of virtual bound states, and 
shown that the element specihc magnetization of Py obtained from XMCD is consistent with 
the ab-initio prediction of VBS formation at Fe sites in the minority band. We expect that 
these states should be the main source of impurity scattering in Py, and can be controlled by 
adjusting the Cn concentration, as was shown by our ab-initio calculations. Theory predicts 
that the density of VBSs at the Fermi level increases with Cn alloying. This should enhance 
the conductivity difference between the majority band and minority band, even though the 
overall reduction in magnetization would certainly counteract this effect. 

We have shown that the exchange interaction between Fe impurities strongly increases 
with Cn alloying. This comes from an increased exchange coupling between localized Fe 
magnetic moments and s — p states at the Fermi level. The VBSs in Py are approximately 
1 eV above the Fermi level. This suggests that electrons excited with an 800 nm wave¬ 
length (~1.5 eV) laser pulse should interact strongly with VBSs in Py, resulting in a large 
exchange interaction between excited electrons and localized Fe moments. As discussed by 
Haag et ah [71H76], the electron-magnon scattering rate depends on the exchange interac¬ 
tion between free electrons and the localized magnetic moments. Hence we surmise that, 
the magnon emission rate should strongly increase at Fe sites for electrons excited ~ 1 eF 
above the Fermi level. Haag et al. m, also showed that the magnon generation is much 
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FIG. 6. Illustration of magnon generation and spin flip processes during demagnetization of Py. 
Free-electrons in the minority band will quickly scatter on VBSs and generate magnons. Electrons 
in the majority band can go to the minority band through a spin-flip process, which is responsible 
for the angular momentum transfer to the lattice. The process can repeat until the electrons reach 
their ground state. 

faster in pure elemental Fe, as compared to Ni, and the magnon emission rate is higher than 
the Elliot-Yafet spin flip rate for both elements. They suggest that, even though magnon 
emission does not change the magnetization due to the concurrent spin-flip process, it pro¬ 
motes s-p electrons from the minority to the majority band, which eventually decreases 
the magnetization via electron-phonon scattering. Similar magnon generation has also been 
compared to pure spin-flip excitations by Schellekens and Koopmans|77]. where they found 
an enhanced demagnetization rate in the former case. This suggests that a qualitatively ac¬ 
curate description of ultrafast demagnetization, must include the effect of electron-magnon 
scattering. 

Enhanced magnon generation at Fe sites should strongly influence the ultrafast demag¬ 
netization processes for Py and Py-Cu alloys. As mentioned in the introduction, Mathias et 
al.flT] observed a signihcant delay between Fe and Ni demagnetization in Pyi_ 2 ,Cu 3 ;. The Cu 
alloying decreases the average exchange interaction between magnetic elements in Py and 
increases the delay of the Ni demagnetization. In Fig. [^we illustrate the proposed demagne¬ 
tization process in Py-Cu alloys. In the minority band, free electrons will rapidly scatter on 
VBSs, and generate magnons (represented by a flipped d-electron), while populating the ma- 







17 


jority band. When this magnon generation is faster than the spin-flip time (r^), there will be 
an accumulation of majority spins. This imbalance between majority and minority spins will 
enhance the demagnetization rate by Elliot-Yafet spin-flip scattering, which is responsible 
for the transfer of angular momentum from the spins to the lattice. Since this process will 
effectively increase the magnetization of the free-electrons, the resulting demagnetization is 
purely a result of magnon generation. Therefore it is critical that the magnon decay time 
{Tmagnon) IS loug euough to accumulate magnons. Commonly, Tmagnon = !/(« ■ 47r/), where 
a ~ 0.01 is the damping parameter and / ~ 1 THz is the frequency region for short wave 
length magnons [75]. This gives Tmagnon ~ 8 ps, which is more than an order of magnitude 
longer than demagnetization times in Py-Cu. Note that the free electrons which transfer the 
angular momentum to the lattice, can repeat this process until they have lost their energy 
and reaches the Fermi level, further enhancing the demagnetization rate of the sample. 

Since the magnons will initially be localized at Fe sites, ultrafast demagnetization should 
initially proceed in the Fe sublattice. It is not within the scope of this paper to provide a 
model of the magnon dynamics. However, we can make the following statement: The group 
velocity of long wave length magnons, Vg ~ 2Dk/h^ is proportional to the average exchange 
interaction, i.e. D oc JAvg, where JAvg is the average pair exchange. Thus the time it takes for 
the magnons, with typical group velocities on the order of 10^ m/s for exchange magnons, 
to distribute throughout the lattice is inversely proportional to the exchange interaction 
that decreases with increased Cu concentration, consistent with the observations of 20-80 fs 
delay between Fe and Ni demagnetizations as found by Mathias et ah HU. If the magnon 
generation rate is higher than the Flliot-Yafet spin-flip rate in pure Fe and Ni, as found by 
Haag et al. m, then the impact of VBSs for the non-elemental resolved demagnetization 
rate should be negligible and the effect of VBSs should only manifest as a demagnetization 
delay between the different elements. 


VIII. CONCLUSIONS 

We have have shown that the decreased magnetization of Py, when alloyed with Cu, is 
the result of Ni d-state localization, in addition to the hybridization of Ni and Cu d-bands. 
Only a small effect from d-band hlling (~6% decrease of Ni 3d-holes with 55% Cu alloying), 
as suggested by the RBM, is observed. Fe in Py and Py-Cu forms VBSs in the minority 
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band above the Fermi level. These states ensure that the Fe majority band remains hlled 
for all Cu concentrations, resulting in a high Fe moment of 2.6 /i^. The intensity of the 
VBS increases with increasing Cu concentration since the average number of host valence 
electrons increases. The increased VBS intensity at the Fermi level increases the scattering of 
minority conduction electrons. Also, this scattering of free electrons on VBS states causes a 
strong exchange coupling between Fe sites, similar to an RKKY interaction. This is observed 
as an increase in the Fe-Fe exchange interaction obtained from ab-initio calculations. In the 
discussion section we propose a mechanism for ultrafast demagnetization in Py and Py- 
Cu, which explains the different responses found for Fe and Ni. The VBSs should have a 
signihcant impact on the probability of generating magnons, which are initially localized at 
Fe sites. The time-delay between the Fe and Ni demagnetization is proportional to the time 
it takes for these magnons to distribute throughout the Py lattice. The magnon generation 
also creates an imbalance in occupation of the spin-dependent s — p bands, which enhances 
the demagnetization rate. 
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X. APPENDIX A: ANDERSON IMPURITY MODEL AND EXCHANGE 

INTERACTION 

Since the VBS formation of Fe is significantly above the ci-band of Ni, it is justihable 
to discuss VBSs using Andersons model, even though the model is only strictly applicable 
for magnetic atoms in a non-magnetic host[HS]. The Anderson impurity model describes 
the mixing between bound states and free electrons, rather than the actual formation of 
the bound states. The half width at half maximum (A) of the VBS, in this model, is a 
direct consequence of the mixing between d— and free electron states, A = n {V^^)ayp{e), 
where Vkd is the mixing parameter and p(e) is the free electron density at the VBS energy 
position. inS] Using A = 0.7 eV and p(e) = 0.046 states /eV /atom, as found from our ab- 
initio calculations for Pyo. 8 Cuo, 2 , we obtain a mixing parameter of Vkd = 2.2 eV. This is in 
the region (2-3 eV) that was predicted by Anderson for 3d-metals. Also, the width remains 
almost constant with Cu alloying, which is consistent with a nearly constant s,p DOS in 
the VBS energy region, found by our ab-initio calculations. 

According to Anderson [65], for magnetic impurities in a non-magnetic metal, the condi¬ 
tion for a magnetic state is where U and J are the effective Coulomb and 

exchange interactions, respectively. We have used U = A eV and J = 0.75 eV which are 
appropriate for Ni in Anderson impurity model [7^ ISO] . The fractional hlling of the d-band 
(uc) is directly given by the number of holes in the non-magnetic state {n,, = 0.66 for Fe 
and 0.87 for Ni). We hnd that Fe clearly satishes the condition for band splitting (12.5 > 8) 
while Ni is far from satisfying it (12.5 ^ 40). Note that A is likely broadened by the CPA 
treatment of the DOS as well as the crystal held splittings, however, even a 3 times smaller 
A will keep the Fe spin split while Ni becomes non-magnetic. 

In the context of Anderson theory it is possible to estimate the exchange interaction 
between magnetic impurities in a metallic host. In this model the exchange interaction be¬ 
tween impurities is mediated by a RKKY-like interaction. Since VBSs are strong scattering 
centers for conduction electrons they can induce a spin polarization of the latter. This po¬ 
larization will mediate the effective pair exchange coupling between the localized moments. 
The exchange between a localized moment and conduction electrons can be approximated by 
'^d-k = , where is the energy of the VBS peak relative to the Fermi level, 

S is the spin of the Fe impurity. However, if one assumes that the scattering electrons (e*,) 
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Vt , Vi (eV) 


FIG. 7. The electron displacement for both minority and majority bands as functions of and 
V-i-, respectively. The electron displacement for different impurities (Cu, Co, Fe, Mn and Cr) in a 
Ni host, is indicated. The corresponding magnetic moments of the impurities are also listed. 

are above the Fermi level and also that there are only VBSs in one of the spin bands, then 
Jd-k = The interaction between Fe moments is then proportional to 


XI. APPENDIX B: CHARGE DISPLACEMENT 

There are essentially two types of VBSs. The hrst corresponds to localized states of 
imparities in a free electron metal, e.g. Fe imparity in A1 host, which is the case described 
by Anderson [HS]. The second case corresponds to changing the energy and behavior of 
localized states that already exists in a host (the 3d band) by adding an imparity, e.g. Fe 
imparity in Ni host, which can be described by the imparity potentials and I 4 in the 
minority and majority bands, respectively. In general, the imparity potential is determined 
by the Coulomb and intra-exchange energies of the impurity atom[62]. The distribution of 
this impurity potential between majority and minority bands is strongly dependent on the 
difference in the number of valence electrons between the host and impurity. In turn, this 
difference will determine if the impurity will be ferromagnetically or anti-ferromagnetically 
aligned to the host. By using the formalism of Campbell and Gomes [62], and Clogston|63]. 
it is possible to use the Ni host DOS and the valence electron difference between the host 
and impurity to self-consistently solve for the impurity potential distribution. In Fig. 
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we present such calculations for various impurities in Ni. We show these calculations solely 
for illustrative purposes. The calculations are identical to what was done by Campbell 
and Gomes[H3) except that a more realistic DOS obtained via ab-initio calculations was 
used. The charge displacement on the left axis corresponds to how much charge has been 
displaced from below the Fermi level to above it due to the impurity potentials Vi and Vf in 
the minority and majority band, respectively. For the majority band, there is a threshold 
before the impurity potential I 4 can displace any electrons. For sufficiently large impurity 
potentials, the majority band can displace more electrons than the minority band. This is 
why Fe and Co, which have a small difference in the number of valence electrons relative 
to the Ni host, will only push minority states above the Fermi level. For larger electron 
displacements, as for Cr, the impurity potential mainly goes to the majority band, which 
makes Cr antiferrromagnetically aligned with the Ni host. (Note that quantitative DFT 
calculations [36] predict a substantial impurity potential in the minority band as well for 
Cr). 


[1] B. Baek, W. H. Rippard, M. R. Pufall, S. P. Benz, S. E. Russek, H. Rogalla, and P. D. 
Dresselhaus, Phys. Rev. Applied 3, 011001 (2015). 

[2] P. Holody, W. C. Chiang, R. Loloee, J. Bass, W. P. Pratt, and P. A. Schroeder, Phys. Rev. 
B 58, 12230 (1998). 

[3] S. Dubois, L. Piraux, J. M. George, K. Ounadjela, J. L. Duvail, and A. Fert, Phys. Rev. B 
60, 477 (1999). 

[4] K. Ando, Y. Kajiwara, S. Takahashi, S. Maekawa, K. Takemoto, M. Takatsu, and E. Saitoh, 
Phys. Rev. B 78, 014413 (2008), 

[5] Y. Niimi, H. Suzuki, Y. Kawanishi, Y. Omori, T. Valet, A. Fert, and Y. Otani, Phys. Rev. B 
89, 054401 (2014), 

[6] B. F. Miao, S. Y. Huang, D. Qu, and C. L. Chien, Phys. Rev. Lett. Ill, 066602 (2013). 

[7] M. Abd El Qader, R. K. Singh, S. N. Galvin, L. Yu, J. M. Rowell, and N. Newman, Applied 
Physics Letters 104, 022602 (2014), 

[8] D. Petrovykh, K. Altmann, H. Hochst, M. Laubscher, S. Maat, G. Mankey, and F. Himpsel, 
Applied Physics Letters 73, 3459 (1998). 













22 


[9] E. Villamor, M. Isasa, L. E. Hueso, and F. Casanova, Phys. Rev. B 88, 184411 (2013). 

[10] M. Haidar and M. Bailleul, Physical Review B 88, 054417 (2013). 

[11] J. Slater, Journal of Applied Physics 8, 385 (1937). 

[12] M. Collins and D. Wheeler, Proceedings of the Physical Society 82, 633 (1963). 

[13] M. Getzlaff, Fundamentals of magnetism (Springer, 2008). 

[14] Y. Kakehashi, Modern Theory of Magnetism in Metals and Alloys (Springer, 2012). 

[15] C. Kittel, Introduction to Solid State Physics (Wiley, New York, 1976). 

[16] R. O’Handley, Modern magnetic materials (John Wiley and Sons, New York, 2000). 

[17] D. Jiles, Introduction to Magnetism and Magnetic materials (CRC press, Boca Raton, 1998). 

[18] A. E. Berkowitz and E. Kneller, Magnetism and Metallurgy (Academic press. New York, 1969). 

[19] D. J. Craik, Magnetism: Principles and Applieations (Wiley, New York, 1995). 

[20] B. D. Cullity and C. D. Graham, Introduction to Magnetic Materials (Wiley, New Jersey, 
2009). 

[21] E. du Tremolet de Lacheisserie, ed.. Magnetism I, Fundamentals (Springer, 1999). 

[22] S. V. Vonsovskii, Magnetism (J. Wiley, 1974). 

[23] R. Skomski and J. Coey, Permanent Magnetism (CRC press, 1999). 

[24] R. M. Bozorth, Ferromagnetism (Wiley-IEEE Press, 1993). 

[25] J. Cable and T. Hicks, Physical Review B 2, 176 (1970). 

[26] P. Dederichs, R. Zeller, H. Akai, and H. Ebert, Journal of Magnetism and Magnetic Materials 
100 , 241 (1991). 

[27] I. Turek, J. Kudrnovsky, V. Drchal, and P. Weinberger, Phys. Rev. B 49 , 3352 (1994). 

[28] R. Jacobs, E. Babic, and J. Xanthakis, Journal of Physics F: Metal Physics 15 , 1941 (1985). 

[29] A. Vernes, H. Ebert, and J. Banhart, Physical Review B 68, 134404 (2003). 

[30] P. James, O. Eriksson, B. Johansson, and I. A. Abrikosov, Phys. Rev. B 59 , 419 (1999). 

[31] T.-S. Choy, J. Chen, and S. Hershfield, Journal of Applied Physics 86, 562 (1999). 

[32] P. E. Mijnarends, S. Sahrakorpi, M. Lindroos, and A. Bansil, Phys. Rev. B 65 , 075106 (2002). 

[33] E. Smirnova, I. Abrikosov, B. Johansson, Y. Vekilov, A. Baranov, V. Stepanyuk, W. Hergert, 
and P. Dederichs, Physical Review B - Condensed Matter and Materials Physics 59 , 14417 
(1999). 

[34] C. Paduani and J. C. Krause, Journal of Applied Physics 86, 578 (1999). 

[35] P. A. Serena and N. Garcia, Phys. Rev. B 50 , 944 (1994). 








23 


[36] R. Zeller, Journal of Physics F: Metal Physics 17 , 2123 (1987). 

[37] V. S. Stepanyuk, R. Zeller, P. H. Dederichs, and I. Mertig, Phys. Rev. B 49 , 5157 (1994), 

[38] P. Yu, X. F. Jin, J. Kudrnovsky, D. S. Wang, and P. Bruno, Phys. Rev. B 77, 054431 (2008). 

[39] N. Y. Moghadam and G. M. Stocks, Phys. Rev. B 71 , 134421 (2005). 

[40] J. Minar, S. Mankovsky, Sipr, D. Benea, and H. Ebert, Journal of Physics: Condensed Matter 
26, 274206 (2014), 

[41] S. Mathias, C. La-O-Vorakiat, P. Grychtol, P. Granitzka, E. Turgut, J. Shaw, R. Adam, 
H. Nembach, M. Siemens, S. Eich, C. Schneider, T. Silva, M. Aeschlimann, M. Murnane, and 
H. Kapteyn, Proceedings of the National Academy of Sciences of the United States of America 
109 , 4792 (2012). 

[42] S. Gunther, C. Spezzani, R. Ciprian, C. Grazioli, B. Ressel, M. Coreno, L. Poletto, P. Miotti, 
M. Sacchi, G. Panaccione, V. Uhlif, E. E. Eullerton, G. De Ninno, and C. H. Back, Phys. 
Rev. B 90 , 180407 (2014), see Supplemental. 

[43] Private communication: Somnath Jana, Joachim Terschluesen, Robert Stefanuik, Johan 
Soderstrom, and Olof Karis (Unpublished, Uppsala University, Sweden). 

[44] A. Eschenlohr, Ph.D. thesis, Universitiit Potsdam (2012). 

[45] B. Koopmans, G. Malinowski, F. Dalla Longa, D. Steiauf, M. Fahnle, T. Roth, M. Cinchetti, 
and M. Aeschlimann, Nat Mater 9 , 259 (2010), 

[46] K. Carva, M. Battiato, D. Legut, and P. M. Oppeneer, Phys. Rev. B 87, 184425 (2013). 

[47] R. Nakajima, J. Stohr, and Y. U. Idzerda, Phys. Rev. B 59, 6421 (1999). 

[48] H. Ebert, D. Kodderitzsch, and J. Minar, Reports on Progress in Physics 74, 096501 (2011). 

[49] http://ebert. cup.uni-muenchen.de/SPRKKR. 

[50] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996), 

[51] C. T. Chen, Y. U. Idzerda, H.-J. Lin, N. V. Smith, G. Meigs, E. Chaban, G. H. Ho, E. Pellegrin, 
and F. Sette, Phys. Rev. Lett. 75, 152 (1995). 

[52] T. N. A. Nguyen, R. Knut, V. Fallahi, S. Chung, Q. T. Le, S. M. Mohseni, O. Karis, S. Pered- 
kov, R. K. Dumas, C. W. Miller, and J. Akerman, Phys. Rev. Applied 2, 044014 (2014), and 
references within. 

[53] M. F. Collins and J. B. Forsyth, Philosophical Magazine 8, 401 (1963), 

[54] R. A. Medina and J. W. Cable, Phys. Rev. B 15, 1539 (1977). 

[55] A. Aldred, B. Rainford, T. Hicks, and J. Kouvel, Physical Review B 7 , 218 (1973), 











24 


[56] J. Tersoff and L. M. Falicov, Phys. Rev. B 25, 4937 (1982). 

[57] J. Vogel and M. Sacchi, Phys. Rev. B 53, 3409 (1996), 

[58] B. Glaubitz, S. Buschhorn, F. Briissing, R. Abrudan, and H. Zabel, Journal of Physics 
Condensed Matter 23 (2011). 

[59] J. Stohr, Journal of Magnetism and Magnetic Materials 200 , 470 (1999), 

[60] G. Meitzner, D. Fischer, and J. Sinfelt, Catalysis Letters 15, 219 (1992). 

[61] N. Stefanou, A. Oswald, R. Zeller, and P. H. Dederichs, Phys. Rev. B 35, 6911 (1987), 

[62] I. Campbell and A. Gomes, Proceedings of the Physical Society 91, 319 (1967), 

[63] A. M. Clogston, Phys. Rev. 125, 439 (1962). 

[64] “Note that the total number of 10 d-states needs to be conserved at the impurity site regardless 
of the number of valence electrons.”. 

[65] P. Anderson, Physical Review 124, 41 (1961). 

[66] J. Friedel, II Nuovo Cimento 7, 287 (1958), 

[67] P. Mohn, Magnetism in the Solid State An Introduction (Springer-Verlag, 2003). 

[68] C. Andersson, B. Sanyal, O. Eriksson, L. Nordstrom, O. Karis, D. Arvanitis, T. Konishi, 
E. Holub-Krappe, and J. H. Dunn, Phys. Rev. Lett. 99 , 177207 (2007), 

[69] J. Stohr and H. C. Siegmann, Magnetism: From Fundamentals to Nanoscale Dynamics 
(Springer, 2006). 

[70] J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966), 

[71] J. Schrieffer, Journal of Applied Physics 38, 1143 (1967), 

[72] K. Yosida, Phys. Rev. 106, 893 (1957). 

[73] C. L. Chien, S. H. Liou, D. Kofalt, W. Yu, T. Egami, and T. R. McGuire, Phys. Rev. B 33, 
3247 (1986), 

[74] M. Haag, C. Illg, and M. Eahnle, Phys. Rev. B 90 , 014417 (2014), 

[75] M. Haag, C. Illg, and M. Eahnle, Physical Review B - Condensed Matter and Materials 
Physics 87 (2013), cited By (since 1996)1. 

[76] C. Illg, M. Haag, and M. Eahnle, Phys. Rev. B 88, 214404 (2013). 

[77] A. J. Schellekens and B. Koopmans, Phys. Rev. Lett. 110, 217204 (2013). 

[78] V. Kambersky and C. E. Patton, Phys. Rev. B 11, 2668 (1975). 

[79] E. §a§ioglu, C. Eriedrich, and S. Bliigel, Phys. Rev. B 83, 121101 (2011). 














25 


[80] B. Surer, M. Troyer, P. Werner, T. O. Wehling, A. M. Lauchli, A. Wilhelm, and A. I. 
Lichtenstein, Phys. Rev. B 85, 085114 (2012), 



